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Abstract. 

We give a general expression for the normally ordered form of a function F[w(a, a')] 
where w is a function of boson creation and annihilation operators satisfying [a, a'] = 1. 
The expectation value of this expression in a coherent state becomes an exact generat- 
ing function of Feynman-type graphs associated with the zero-dimensional Quantum 
Field Theory defined by F(w). This enables one to enumerate explicitly the graphs 
of given order in the realm of combinatorially defined sequences. We give several 
examples of the use of this technique, including the applications to Kerr-type and 
superfluidity- type hamiltonians. 
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In the normally ordered form of a function F(a, a)) of boson creation and 
annihilation operators all the annihilation operators are moved to the right using the 

= 1. The importance of the normal form, denoted by 



a, a 



commutation relation 
J\f[F(a, a')] and satisfying F(a, a*) = Af[F(a, a')], is evident, as with it the expectation 
values of F(a, a') can be easily evaluated in such canonical states as the vacuum |0) and 
coherent states \z) = e~' z ' 2 / 2 J^Lo^Vv^iF l n )> ( z com pl ex an d a)a\n) = n\n)). The role 
of the normal form in Quantum Field Theory (QFT) is pre-eminent [1, 2]. 

In this work we consider functions F(w) that involve operators w(a, at) in the 
form of a product of positive powers of a) and a, and powers of (a + a*), although the 
formulae derived below are valid for more general i&'s. Our considerations are based on 
the following operational property of formal power series. Let f(x) = Xl^Lo /n 2 -"/ 71 ' an d 
d( x ) = Y^=o9nX n /n\ be two formal power series, also called the exponential generating 
functions (egf) of sequences {f n }^=o and {g n }'^ ) =0 , respectively. Then it can be verified 
that 



' "" dx 



9 



x=0 



x=Q 



oo 

E 

n=0 



fn 



A™ 
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n\ 



which we shall call the product formula (PF). This implies the following property 
satisfied by F(Xw) with indeterminate A: 



F(Xw) 



A^ 

dx 



x=0 



which by taking the normal form of both sides becomes 

d \ 



J\f[F(Xw) 



A-^ 1 M{e a 
dx J 



(1) 



(2) 



x=0 



Note that in Eq.(2) a separation has been achieved between the functional aspect 
(defined by F) and the operator aspect (defined by w) of the normal ordering. 
Conventionally we implement M{e xw ) by using the auxiliary symbol : : with M{e xw ) = : 
G$(x, a, a') : , where under the symbol : : the function G^(x, a, a') is normally ordered 
assuming that a) and a commute [3, 4]. Then Eq.(2) becomes 

M[F{\w)] = F : Gu,(x,a,a*) : 

Note that the expression of Eq.(3) arises in the evaluation of the partition function Zp 
for the system defined by the hamiltonian 7i(w) 

1 



(3) 



Tr e 



-m{w) 



TV 



d l z 



x=0J 



taking the trace over the coherent state representation, (3 = (/c B T) _1 [5]. The problem 
of finding Af[F(\w)] reduces to that of finding Af(e xw ), still however a non-trivial 
task, vide the classical references [3, 4, 6]. We have recently found expressions for 
Gw(x, a, at) for several types of operators w of the form W( rjS ) = (o)) r a s [7] as well as for 



w 



(r,s) 



n 



M 



k=l W (r k ,s k ) 



with r, s, Tfc, Sk positive integers [8]. 



At this point it is already possible to relate Eq.(3) to enumerative formulae for 
Feynman-like graphs in QFT [9]. Assume that our formal power series F(x) can be 
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written in the form F(x) = exp (Ylm=i ^m^r), and we similarly assume that we may 
define operators Vn W \a,a)) by 



: a, a*) := : exp ( ^vfW) — 



,n=l 



TV. 



(4) 



Explicit examples [7] from which the operators Vn (a, a') may be read off include: 

w = a) a, A/" [exp (rca^a)] = : exp [ata(e x — 1)] : , (5) 

w = (a)) r a, ATfexp (x(a^ r a)} = 
= : exp 



(6) 



with r = 2, 3 . . . , as well as more involved expressions for other w(a, a'). Thus, Eq.(3) 
may be written as 

J\f[F(\w)] = exp (E~=i ^A-^) • : exp fe~ , ) = n • 

\ / x=0 

We eliminate the operators a and at by taking the matrix element of Eq.(7) in the 
coherent state \z) and using a\z) = z\z). This yields 

{z\N[F{\w)]\z) = exp (£~ =1 ^\ m £.) ■ exp (E^^(M*)j 

By specifying z = 1 in Eq.(8), defining \4 ( ^(1,1) = V^, V = {K W }^ and 
L = {L m }^ =1 we obtain 

Z(L,V,X) = (l\Af[F(Xw)]\l) 



(7 



x=0 



exp ^ 



m! rfx m / \ ^ n n\ 

m=l / \n=l 



(9) 



a;=0 



which is essentially the counting formula cited by Bender et al. [9] . Due to the symmetry 
of the PF , we have Z(L, V, A) = Z(V, L, A), which may facilitate the calculations. 
Furthermore it can be demonstrated that for all the forms of w used here the sequence 
V consists of positive integers. The formula Eq.(9) was employed in [9] as an enumerative 
tool for counting the Feynman-like graphs in zero- dimensional QFT models, where the 
values of all Feynman integrals are equal to one. Our derivation sheds light on its 
quantum origin by tracing back its sources to the boson normal ordering problem. By 
specifying the sets L and V one can attempt to produce a (in general divergent) power 
series expansion in A: 

Z(L,V,A) = ^A n (L,V)- (10) 



n=0 



in which A n (L, V) can be related to known objects. To see that, recall the definition 
of the multivariate Bell polynomials M(h,u) related to a function h(x) = X^L^n^r 

through (h = {Mn=l) 



uh[x) 



= E ^ E ukB **<v = E ^ B «( h > «)» 



m 

n=0 k=l 



(11) 



n=0 
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where the coefficients of the expansion M n (h,u) = ^^ =1 -u fc B nfc (h) depend only on 
hi, ... , h n . We refer to [10, 11] for further properties of B n (h, u). 

With B n (f) = B n (f, 1) we see that the coefficients A n = A n (L, V) factorize 

A n = B n (L) ■ B n (V) (12) 

which for given L and V can be worked out (see below). 

The utility of Eqs.(9) and (12) goes beyond the specific definition of initial w, and 
this is the philosophy of Ref. [9] where it was suggested that L and V could be treated 
as initial input for QFT models. From this perspective Eqs.(9) and (12) provide the 
starting point for a Feynman-like graph representation of the coefficients A n in Eq.(10), 
where A n counts the number of graphs with n labelled lines. The graph construction 
rules are: a line starts from a white dot, the origin, and ends at a black dot, the vertex. 
We further associate strengths \4 with each vertex receiving k lines and multipliers 
L m with a white dot which is the origin of m lines. Counting such graphs consists in 
calculating their multiplicity due to the labelling of lines and the factors L m and 

We now specify L and V and give some examples of the explicit evaluation of A n 
along with the explicit graph representation: 

Example 1: L\ = 1, Lm = 1 (M > 1), and L rn = otherwise, giving the function 
F(x) = exp(x + x M I Ml); Vn W ^ = 1 for n — 1, 2 . . ., which arises from the string w = a^a, 
see Eq.(5). This corresponds to the normal ordering problem J\f[exp(\a^a + j^{a)a) M )\. 
Note that the case M = 2 describes the normal ordering of the exponential of the Kerr- 
type hamiltonian [12] H = \a^a{l + ^a)a). Using the definition of the two variable 
Hermite-Kampe de Feriet polynomials H n M \x,y) (see [13] and references therein) 

Y, Hi n M) ^y)h = eXt+ytM > ( 13 ) 



n=0 



y.n — Mr „.r 



where H n \x,y) = n\ ^jZ ^ n _ Mr y rV F(x) can be expanded as 

1 x n 
M! ra! 



F(x) = e^ = J2H n -\l,±f-. (14) 

n=0 



Eq.(10) yields A n = H n MS> (l, -^y) ■ B n , where the Bell numbers B n are defined 
through their egf: exp(e x - 1) = Er?=o 5 ™S t 9 ' 10 ' Observe that for M = 2, 

H [ n\l,\) = (75)" H n (-j=^ = 1,2,4,10,26,76,232,... are the involution numbers 
[10] expressible using Hermite polynomials H n (x). The initial terms of A n for M = 2 
are: 1,4,20,150,1352,15428,..., see Fig.(l), and for M = 3: 1,2,10,75,527,6293,..., 



etc. Note that whereas B n counts all the partitions of an n-set, Hn M \l, -gj) counts 
partitions of an n-set into singletons and M-tons. 

Example 2: L m = m for m = 1,2... , giving rise to F(x) = exp(E^ =1 m^f) = 

exp(xe x ) = E^Lo^™^"' wn ere /„ = EfcLo ( ^ ) k n ~ k are idempotent numbers 

[10]. Again choosing = 1, n = 1,2..., with w = a^a, gives A n = 
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n=1: 



1 = 1 



> < 



n=2: OCT Q 



1 + 1 + 1 + 1=4 



n=3: 




1 + 3+1 



() • 



< 




+ 3 +6 + 3= 20 



Figure 1. Lowest order Feynman-type graphs for Example 1 with n — 1,2,3 lines. 
The number below each graph is its multiplicity. 



I n ■ B n = 1,6,50,615,10192,214571... . This corresponds to normally ordering 
AT[exp(A(a t a)e A ( ata ))]. 

Example 3: L\ = 0, L m = 1 for m = 2, 3 . . . , leading to F(x) = exp(e x — 1 — x) = 
Y^=o^n j, where B$p are restricted Bell numbers which are defined as counting 
partitions without singletons. (Note that Bn = - YlkLo )■ Here we choose = 
n\, n = 1,2... , derived from the string w = (a^) 2 a, and producing via exp (^z - ) — 



£~=o 5 « f , ( see E( l-( 37 ) of [7]), 4» = M iJ • B\T ' = 0, 3, 13, 292, 5511, 166091 . . , see 
Fig. (2). This corresponds to the normal ordering of exp ^e x( - at ^ a — 1 — A(a^) 2 aj. 

Example 4: L 2m = 2m, L 2m +i = for m = 0, 1, 2 . . . , giving F(x) = exp(x sinh(x)). 
If, as in Examples 1 and 2, = 1, n = 1,2..., w = a^a, then by defining 

the idempotent polynomials I n (t) = J2k=o (Ju) k n ~ k t k we obtain A n = 1$ ■ B n = 

0,4,0,240,0,49938,0,24608160,0... , where iP = ^ =0 ( n ) (-l)fc/ fe (-I)/ n _ fe (I), 

yielding J\f[exp(Xa^a sinh(Aa'a))]. 

Example 5: In the last example we shall treat the function w = a + a\ using 
F(x) = e xM ' M \ M = 1,2,3,.... First observe that N{e™) = : G^{x,a,a^) : = 
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n=2: 



<: 



1 * 1 + 7*2 



1x1 + 3x2 + 1 x 6 



13 



n=4: 





3*7 + 1 x 1 



1 

2 
I 

72 x 2 



<: 



O * 2 



6*2 




6*2 



O • ? 



3*4 + 6*4 



3*4 




72 * 6 + 4 * 6 + 3 * 24 + 1 * 24 = 292 



Figure 2. Lowest order Feynman-type graphs for Example 3 with n — 2, 3,4 lines. 
The number below each graph is (multiplicity) x W^ivertex factor Vk = k\). Numbers 
at black dots (vertices) are vertex factors. 



: e x2 / 2 e x ( a + at ) : which is a consequence of the Heisenberg algebra. It follows that 
Vi(a,al) = a + a\ (a, at) = 1 and Vn\a,at) = for n > 2, see Eq.(4), giving 
V = {2, 1, 0, 0, ...} and L = {5 m ,M}m=i- Let us define the modified Hermite polynomials 
K{x) = {-^ n H n [^j and then exp (2x + x 2 /2) = £~ Q ^x". Usin S Eqs.(10) 
and (12) we get 

( \ M ol M \ ( x 2 \ 
Z M (L, V, A) = ex P {jw.lxMj ' 6XP { 2X+ YJ _ 

-^[m ■ (15) 

n=0 v 7 

Starting with the simplest case M— 1, the function Zi(L, V, A) = exp (2A + A 2 /2) gives 
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n=2: 



<: 



1 * 4 



1 * 1 



n=4: 



<: 
<: 



3 * 16 



6 * 4 



<: 



12* 4 



+ 3*1 



6 * 1 



Figure 3. Lowest order Feynman-type graphs for Example 5 with n = 2,4 lines. The 
number below each graph is (multiplicity) x Y\ k {vertex factor Vk = 2<5fc.i + 5k,2)- 



A n = h n {2) = 1,2,5,14,43,142,499,1850,..., n = 0,1,2,... . The series of Eq.(15) 
can also be written down in closed form for M — 2, corresponding to a single mode 
superfluidity-type hamiltonian TC ~ (a + a)) 2 [14], and for M = 3 [15]: 

7 (J V U V^^( 2 ) f X2 Y 1 ( \ 

ML,v,\)=j:—(^) =jT^W^ exp [T^)> (16) 

n=0 v 7 

(l-^A 3 ) 1 ^ 2r °V6'6' '2(1-0A 3 )V' 1 J 

where 0(A) = 1 ~ v/ ^ 3 " 4A3 and 2-^0 is the generalized hypergeometric function (the Airy 
function, [16]). In these examples Z\ and Z2 are convergent series in A while Z% is 
not - summation of this series is understood in the generalized sense. From Eq.(15) 
we can read off the values of A n : Amu = (my^l^ Mni^) and zero otherwise, giving for 
M = 2, A 2n = 1, 5, 129, 7485, 755265, 116338005,'..., see Fig.(3). Note, that whenever Z 
is known in closed form the equation Z(\) = exp (a^) Z[x)\ x _ q leads immediately to a 
set of graphs for which L m = <5 mj i. Thus for M = 2, with Eq.(16) we have the following 
alternative descriptions: a) L m = 5 m)2 ; V\ = 2, V% — 1, V n> 2 = and b) L m = S m> i; 
Vm = (4n + l)(2n — 1)! . However even if Z is not known explicitly method a) leads to 
a simple, alternative, graphical description using Eq.(15). 

In conclusion we see that the technique described herein and hinging on Eq.(9) 
leads to a combinatorial and graphical description of many physical systems. 
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